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I.  INTRODUCTION 

Not  only  is  the  physical  interpretation  of  the  radiation  of  the  acoustic  waves  in 
more  than  one  dimension  complex,  but  the  effects  of  varying  frequency  and  shading 
are  hard  to  visualize  when  considering  directionality  and  frequency  properties.  As 
Wood  implied  in  his  book,  A  Textbook  of  Sound  [Ref  1:  p.449],  the  directional  factor 
is  very  important  in  many  acoustic  applications  such  as  speech,  hearing  and  array 
design. 

The  simplest  array  one  can  consider  is  the  beam  pattern  of  a  dipole  radiating  in  two 
dimensions.  References  to  the  case  of  the  simple  dipole  can  be  found  in  virtually  every 
textbook  in  acoustics.  For  a  fairly  recent  example,  see  Malecki's  book  [Ref  2:  p.  199]. 
Currently,  the  mathematical  modeling  of  the  simple  acoustic  arrays  has  been  enhanced 
with  the  advent  of  the  sophisticated  computer  software.  Specifically,  Ziomek 
[Ref  3:  p. 138]  and,  Embleton  and  Thiessen  [Ref  4:  p. 1124]  have  investigated  N  point 
source  arrays,  considering  both  amplitude  and  phase  weighting.  Although  the 
mathematics  and  computer  hardware  exist,  the  visualization  of  the  three-dimensional 
aspects  of  these  arrays  is  difficult.  By  use  of  the  computer  graphics  software  the  ability 
to  interpret  the  significance  of  shading  and  phasing  can  be  fully  developed  [Refs.  5,6]. 

To  date,  visualization  of  the  three-dimensional  beam  pattern  has  been  needed  to 
understand  the  directional  factor.  Numerical  computer  output  was  essential  to  draw 
three-dimensional  graphs.  By  changing  two  parameters  (in  three-dimensions,  those 
could  be  X  and  y  or  0  and  (p,  etc.),  different  generated  values  of  the  beam  pattern,  F, 
are  obtained. 

One  computer  program  to  visualize  the  beam  pattern  in  three-dimensions  (P3D)  and 
four  computer  programs  to  read  values  in  detail  (2DICIR,  2DPCIR,  2DIFIX, 
2DPFIX)  were  written.  These  computer  models  predict  the  directional  factor  and 
overall  beam  pattern  with  or  without  differences  in  the  source  strengths,  inter-source 
distances,  and  phases  of  the  elements  of  the  arrays. 

Typical  geometry  for  the  beam  pattern  prediction  is  shown  in  Figure  2.3  on  page 
16. 

The  following  defmitions  apply  and  will  be  used  throughout: 

1  =  source  strength  of  the  first  source  (source  1) 


A  =  source  strength  of  source  2 

B  =  source  strength  of  source  3 

C  =  source  strength  of  source  4 

Dj  =  distance  from  the  geometrical  center  to  source  1 

D2  =  distance  from  the  geometrical  center  to  source  2 

D3  =  distance  from  the  geometrical  center  to  source  3 

D^  =  distance  from  the  geometrical  center  to  source  4 

(^2  -  phase  difference  between  source  2  and  source  1 

(Pj  =  phase  difference  between  source  3  and  source  1 

(p^  =  phase  difference  between  source  4  and  source  1 

9  =  angle  measured  from  the  z  axis  toward  the  x-y  plane 

<J>  =  angle  measured  counterclockwise  in  the  x-y  plane  starting  from  the  x  axis 

k  =  wave  number 

w  =  angular  frequency 

R  =  distance  from  the  geometrical  center  to  a  certain  point  in  space 

H  =  directional  factor 

F  =  beam  pattern  =  201og^Q|H|  or  lOlog^QH^  expressed  in  dB 

All  distances  are  measured  from  the  origin  of  the  coordinate  system  to  the  sources 
along  the  x  and  y  axes.  Phase  differences  are  based  upon  the  phase  of  the  first  source, 
which  is  assumed  to  be  zero. 

The  goal  of  this  thesis  is  to  develop  a  convenient  methodology  for  the 
investigation  of  the  functional  dependency  of  the  directional  factor  on  the  source 
amplitudes  (A,  B  and  C),  distances  from  the  origin  (Dp  D2,  D3  and  D^),  phase 
differences  ((p^.  ^3  and  (p^),  and  0  and  <I>  as  demonstrated  by  the  five  programs:  P3D, 
2DICIR,  2DPCIR,  2DIFIX,  2DPFIX. 
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II.  THEORY 

A  spatial  configuration  of  simple  sources,  either  discrete  or  distributed,  each  with 
its  own  complex  source  strength,  can  be  used  to  represent  the  quadruplet  source  case. 
The  pressure  at  a  field  point  is  the  sum  of  pressures  produced  by  the  individual  sources 
[Ref!  7:  p.  169].  Before  we  develop  the  quadruplet  source  case,  the  acoustic  doublet 
case  needs  to  be  described  to  provide  a  background  for  the  more  complicated 
quadruplet  case. 

A.       DOUBLET  SOURCES 

A  doublet  source  consists  of  two  point  sources  of  strengths  of  1  and  A,  separated 
from  the  origin,  by  distances  D^  and  D2  and  vibrating  at  the  same  frequency  but  with 
a  phase  difference  (p. 


1) 

f 

V" 

A' 

0                       / 

Figure  2.1     Geometry  of  doublet. 
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From  the  geometry  shown  in  Figure  2.1,  we  have 

cose  =  8j/Dp  (2.1) 

sine  =  e^/Dp  (2.2) 


r^  =  V(R-ei)^  +  ^{ 


=  ^{R  -  D^sine)^  +  (Djcose)^,  (2.3) 

-i-  =  V(l  -  D^sine/R)^  +•  (6jCose/R)^.  (2.4) 

sine  =  C2/D2,  (2.5) 

cose  =  S^/D^,  (2.6) 


Tj  =  V(R  +  £2)^  +  ^2^ 


=  V(R  +  D^sinO)-  +  (D^cose)^  (2.7) 

i-  =  V(l  +  D^sine/R)^  +  (S^cose/R)^.  (2.8) 


The  pressure  at  point  (R,e)  due  to  source  1  is 
P    =  J_  eJ(^  -  ^i) 


'     h 


(2.9) 


where  w  is  the  angular  frequency.  That  due  to  source  2  is 


p    =  ^eK^-^2-^).  (2.10) 

The  total  acoustic  pressure  at  point  (R,e)  is  then  the  sum  of  P,  and  P^, 


p  =  -L  eK^vt  -  kR)(  J.  ,ik(R  -  rj)  +  :1L  .-Mr^  -  R)  +  <Pl  }  (2.11) 

where  R  is  the  distance  from  the  field  point  to  the  origin.   From  (2.4)  and  (2.8),  we  get 

k(R-rp  =  kR(l  -  r^/R) 
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=  kR[l  -  V(l  -  D^sinG/R)^  +  (D^cosG/R)^  ],  (2.12) 

k(r2  -  R)  =  kRCr^/R  -  1) 


=  kR[V(D2sine/R  +  1)^  +  (D2cos0/R)2  -  1  ].  (2.13) 

Therefore,  equation  (2.11)  can  be  rewritten  as 

p  =  -L  eJ(wt  -  kR)f  _5i  gjkR[l  -  ^(1  -  DjSin0/R)2  +  (D^cos0/R)2  ] 
R  r, 

+  _^e-J[^^f>/(D2S^^/R  +  1)^  +  (D2cos0/R)2  -  1  ]  +  cp]  I         (2  14) 

In  the  most  frequently  encountered  cases,  the  observation  of  the  pressure  is  made 
at  distances  greater  than  the  separation  of  the  sources.  Therefore,  it  will  be  useful  to 
derive  the  far  field  (R>  >  (Dj  +  D2))  pressure  equation  from  the  above  equation  (2.14). 
From  the  geometry  shown  in  Figure  2.2,  we  have 

sin0  =  e^/D^  =  £2/02.  (2.15) 

Therefore,  (R  -  r^  and  (r2  -  R)  can  be  replaced  with  DjSin0  and  D2sin0.  Also  we  have 

R  R 


^1  h 


-  1  (2.16) 


so  that  the  pressure  equation  becomes 


p  =  i_  eJ(wt  -  kR)jgjkDjSin0  +  ^g-j(kD2sin0  +  (p)j  ^2.17) 


|P|  =  P(R)H(0)  (2.18) 

where  P(R)  =  (14-A)/R. 
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H(e)  = 


1 


(l  +  A) 


JCkDiSine)  +  ^g-JCkDjSine  +  (p) 


(2.19) 


This  directivity  factor,  H(0)  can  be  written  as 


H(0)  = 


1 


(l  +  A) 


•{[cos(kDjSine)  +  Acos(kD2sine  +  (p)]^ 


+  [sin(kDjSine)  -  Asin(kD2sine  +  (p)n 


l^Yz 


(2.20) 


Figure  2.2    Geometry  of  doublet  (far  field). 
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B.       QUADRUPLET  SOURCES 

Similarly,  a  quadruplet  source  case  can  be  developed  with  four  sources  of 
strengths  1,  A,  B  and  C,  separated  from  the  origin,  by  distances  D^  D^,  D3  and  D^ 
and  vibrating  at  the  same  frequency  with  phase  differences  cp^,  <P3  and  cp^  (Figure  2.3). 


Figure  2.3     Geometry  of  quadruplet. 
From  the  geometry  shown  in  Figure  2,3,  we  have 


X  =  RsinGcosO, 
y  =  Rsin6sin<I), 
z  =  Rcos6, 


R  =  Vx^  +  y^  +  Z' 


2  ,  ,2 


(2.21) 
(2.22) 
(2.23) 

(2.24) 
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where  x,  y,  z  are  rectangular  coordinates  and  R  is  the  distance  from  the  field  point  to 
the  origin. 
And  then 


Tj  =  V(y  +  Dj)^  +  X-  +  Z-.  (2.25) 

By  Taylor's  approximation 

r^  -  R  +  D^y/R  +  Dj2/2R.  (2.26) 

r^  -  R   =  D^y/R  =  D^sinGsinCp,  (2.27) 


r^  =  V(y  -  D^)^  +  x^  +  z-  -  R  -  D^y/R.  (2.28) 

R  -  2  =  D^y/R  =  D^sinesinO  (2.29) 

where  R  is  much  greater  than  D^  and  D2. 
Using  the  same  approach 


r,  =  V(x  -  D,)^  +  y-  +  z-  =  7R^  .  2D.X  +  D,^ 


=  RVl  -  (2D3X/R^  -  Dj^/R^)  .  (2.30) 

rj  -  R[l  -  (1/2)(2D3X/r2  -  Dj^/R^)]  =  R  -  D3X/R  +  D32/2R,  (2.31) 

R  -  r3  =  D3X/R  =  D3sin0cosa)  (2.32) 


1*4  =  V(x  +  D^)^  +  y^  +  z-  -  R  +  D^  x/R,  (2.33) 

r_^  -  R  =  D^sinGcosO  (2.34) 

where  R  is  much  greater  than  D3,  and  D^. 
Therefore,  the  pressure  at  (R,  0,  (p)  is 

p  =  JL  eJ(^vt  -  kR  -  kDjSinGsinO)  +     ^   gj(wt  -  kR  +  kD2sin9sin<I)  -  (^2) 
R  R 

^     ^     gj(wt  -  kR  +  kDjSinGcosO  -  (Pj) 
R 

+  ej(wt  -  kR  -  kD^sinGcosO  -  cp^)  /2  25) 

R  ■  ^     ^ 
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The  directivity  factor, 

H(e  (D)  =  ! I  e-i(^DiSin^sin(D)  +  ^gj(kD2sinesin<D  -  (p-^) 

(1  +  A  +  B  +  C) 

+  ggjCkDjSinGcosO  -  93)  +  (^g-KkD^sinGcosO  +  cp^  j^  ^2.36) 

Expanding  the  complex  tenns  yields 

H(0,<P)  =  {[cos(kD,sin0sin<I>)  +  AcosCkD.sinGsinO  -  9J 

+  Bcos(kD3sin0cos<I)  -  93)  +  Ccos(kD4sin0cos<I)  +  (p^)f 

+  [sin(kDjSin0sin<I))  -  Asin(kD2sin0sin<P  -  (p^) 

-  Bsin(kD3sin0cos<D  -  93)  +  Csin(kD^sin0cos<D  +  (p^)]^}'"'.        (2.37) 

This  is  the  expression  used  for  the  computer  program. 

C.       GRAPHICAL  PRESENTATION  OF  THE  DIRECTIVITY 

1.  T>vo-Dimensional  Rectangular  and  Polar  Plots 

From  equation  (2.37),  H  (0,<I>)  versus  0,  and  H  (0,<I))  versus  O  can  be  plotted 
if  either  O  or  0  respectively-  is  given  constant  value.  In  a  rectangular  plot,  the 
directivity,  F,  which  is  equal  to  201og^Q|H  (0,0)  |  is  plotted  on  one  axis  and  0  or  O  on 
the  other.  An  example  of  this  is  seen  in  Figure  3.3.  In  a  polar  plot,  the  directional 
factor  or  beam  pattern  is  a  magnitude  of  the  radial  distance  from  the  center  of  the  plot 
while  the  appropriate  angle  is  measured  clockwise  around  the  origin.  An  example  of 
this  is  seen  in  Figure  3.4. 

2.  Three-Dlmensional  Cartesian  Plot 

From  equation  (2.37),  H  can  be  calculated  and  plotted  as  a  function  of  both  0 
and  O.  This  provides  the  three-dimensional  surface  with  its  elevation  equal  to 
201og^Q|H  (0,O)|  for  each  0-<I>  pair.  An  example  of  this  is  Figure  3.7.  This  form  of 
presentation  is  difficult  to  interpret  but  can  be  a  very  useful  analytical  tool.  Further 
explanations  are  given  in  the  next  chapter. 

3.  Contour  Plot 

A  contour  plot  gives  lines  of  constant  H  (0,O).  As  an  example,  let  kD^  = 
kD2  =  kDj  =  kD^  =  kD,  (p^  =  0,  (pj  =  (p^  =  tt  and  A  =  B  =  C  =  1.  Equation 
(2.37)  then  can  be  simplified  to 

17 


H(e,a))  =  |cos(kDsin0sin<I))  -  cos(kDsinecosa))|.  (2.38) 


By  setting  H(8,<I))  =   Q,  the  desired  value,  and  sinGsin^  =  Cp  and  sinGcosO  =  C2 

where  [c^i  -^  IC2I,  we  get 

lcos(kDc^)  -  cosCkDc^)!  =  2Q.  (2.39) 

This  transcendental  equation  is  then  solved  for  Cj  and  c->.  Now,  with  c^  and  c^ 
determined,  we  need  to  solve  the  second  transcendental  equation 

tanO  =  -^  (2.40) 

which  was  derived  by  dividing  Cj  by  c,,  for  <[>  and  then,  for  G  using  the  G-O 
relationship.  However,  many  of  c,-C2  pairs  are  possible  solutions  to  the  first 
transcendental  equation.  Therefore,  equation  (2.40)  must  be  solved  for  all  the  different 
Cj-Cj  pairs  by  repeating  the  above  procedure.  This  would  result  in  a  contour  plot  such 
as  Figure  3.10.  This  analytical  approach  was  not  pursued  further  in  this  research 
because  other  algorithms  for  extracting  constant  values  of  F  for  varying  G  and  O  had 
already  been  obtained  and  successfully  applied.  However,  the  utility  of  the  above 
transcendental  equation  approach  may  well  bear  further  investigation,  and  we  will 
pursue  a  little  case-by-case  analysis  for  some  cases  for  which  the  solution  of  equations 
(2.39)  and  (2.40)  is  relatively  straightforward. 
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III.  PROGRAM  DEVELOPMENT 

In  analyzing  the  directional  and  frequency  properties  of  shaded  and  phased 
simple  arrays  several  types  of  display  methods  were  needed.  At  the  first  stage  of 
developing  the  computer  program,  a  HEWLETT-PACKARD  86B  computer  and  a 
HP7090A  MEASUREMENT  PLOTTING  SYSTEM  were  used.  By  using  this 
equipment  we  were  able  to  get  two-dimensional  plots  as  shown  in  Figures  3.1  and  3.2. 
However,  it  took  much  time  in  executing  the  program  and  also  we  could  not  do  three- 
dimensional  plottings. 
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Figure  3.1     2-D  rectangular  plot  (HP  86B). 
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PHI    VS   DIRECTIVITY 


180(DEG)     - 


K*D   =  2*PI 

A=l        B=l        C=l 

P1=0     P2=0     P3=0 


O(DEG) 


Figure  3.2    2-D  circular  plot  (HP  86B). 

To  get  a  complete  understanding  of  the  beam  pattern,  three-dimensional  displays 
viewed  from  various  angles  were  essential.  Consequently,  alternative  methods  of 
display  were  programmed  in  FORTRAN  and  executed  on  the  IBM  3033.  There  are 
several  reasons  why  this  proved  to  be  the  most  desirable  method.  Primarily,  a  three- 
dimensional  graphics  software  package  (DISSPLA)  is  available  for  use  on  the  IBM 
3033.  Secondly,  FORTRAN  generated  the  numerical  output  with  greatly  increased 
efficiency  when  using  the  EXTENDED  FORTRAN  COMPILER  [Ref  8:  p.l8]. 
Finally,  the  IBM  3033  has  many  more  output  devices  such  as  VERSTEC,  TEK618  and 
SHERPA  to  name  a  few.  Therefore,  the  programs  used  in  the  HEWLETT- 
PACKARD  86B  were  rewritten  in  FORTRAN.  The  translation  from  BASIC  to 
FORTRAN  required  some  changes  in  variable  names  along  with  the  use  of  some 
convenient  function  calls.  For  reference,  some  of  the  original  BASIC  programs  are 
attached  in  Appendix  A. 
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Figure  3.3    2-D  rectangular  plot  (DISSPLA). 

In  the  trial  of  analyzing  only  with  numerical  outputs,  BASIC  had  become 
inefiective  in  recognizing  details  due  to  the  very  large  number  of  computed  values.  At 
this  stage  the  graphic  package  DISSPLA,  available  on  the  IBM  3033  computer,  was 
used  in  conjunction  with  these  programs.  The  DISSPLA  package  is  a  library  of 
FORTRAN  subroutines  that  facilitate  data  plotting.  It  is  also  device  independent, 
which  means  that  it  does  not  rely  upon  features  particular  to  any  type  of  graphical 
device  and  has  its  own  extensive  symbol  and  character  generation  routines.  By  being 
device  independent,  programs  of  the  HEWLETT-PACKARD  86B  could  be  translated 
without  difficulty.  Therefore,  DISSPLA  was  used  to  produce  both  two-dimensional 
and  three-dimensional  graphs.  The  data  used  to  compute  graphs  could  be  generated 
within  the  FORTRAN  program  or  read  from  a  data  file.  The  former  method  was  used 
here  because  data  files  took  a  great  amount  of  file  space  due  to  the  large  volume  of 
values. 
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PHI  VS.  DIRECnVITY 


THEn'A=35  DEGREE 
K=10.5 

Dl=0.5  D2=0.5 
D3=0.5  D4=0.5 
P2=0   P3=0   P4=0 
A=l.  B=l.  C=l. 


Figure  3.4    2-D  circular  plot  (DISSPLA). 

Two-dimensional  graphs  were  initially  tried.  By  holding  one  parameter  constant 
and  -with  the  use  of  some  array  variables,  plots  of  directivity  versus  angle  were 
obtained.  Examples  of  two  kinds  of  displays  are  shown  in  Figures  3.3  and  3.4. 
Programs  are  attached  in  Appendix  B. 

To  predict  the  real  shape  of  a  directivity  pattern  in  three-dimensions,  three- 
dimensional  graphics  had  to  be  developed.  In  trying  to  get  some  satisfying  three- 
dimensional  plotting,  all  methods  in  DISSPLA  were  tried.  They  are  SURFUN,  and 
SURMAT  for  regular  matrix,  scattered  points,  and  three-dimensional  vector  dra^v^ing. 

SURFUN  is  a  subroutine  that  draws  the  surface  from  points  in  the  z  direction  in 
conjunction  with  a  grid  in  the  x-y  plane  on  the  condition  that  z=f^x,y)  is  given.  One 
constraint  that  was  faced  in  applying  this  function  was  that  there  was  more  than  one  z 
value  for  a  specific  x-y  coordinate. 
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Figure  3.5     3-D  plot(SURMAT  for  scattered  points). 

The  next  trial  was  SURMAT  for  a  regular  matrix  (for  example,  a  rectangular 
three-dimensional  lattice  of  points).  For  the  same  reason  as  in  applying  the  SURFUN 
subroutine,  this  also  turned  out  to  be  non-applicable  even  though  all  data  points  were 
translated  into  a  regular  matrix. 

The  next  option  tried  was  SURMAT  for  scattered  points  (for  example,  a  three- 
dimensional  irregular  cluster  of  points).  As  seen  in  Figure  3.5  it  also  was  unsuccessful 
due  to  more  than  one  point  on  one  specific  coordinate  of  the  x-y  plane.  The  program 
is  hsted  in  Appendix  C  for  reference. 

The  last  option,  the  three-dimensional  vector  drawing  method  gave  us  a  more 
understandable  display  as  shown  in  Figure  3.6.  However,  it  did  not  show  the  beam 
pattern  clearly  because  all  the  lines  were  visible,  even  from  different  view  angles.  At 
this  stage  the  problem  was  the  removal  of  hidden  lines.  Unfortunately,  that  was  not 
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possible  because  the  starting  and  ending  points  of  the  hidden  lines  could  not  found. 
The  program  is  listed  in  Appendix  D. 


DIRECTIVITY  PATTERN 

<^ 

^ ^^^^^"^^ 

K=10.55 
Dl=0.5  D2=0.5 
D3=0.5  D4=0.5 
VIEW  ANGLE=  45,  20 
P2=0  P3=180P4=180 
A=L  B=l.  C=l. 

^^K 

^ic^^»^^^^'^wv^« 

^^^E 

L=^ 

Figure  3.6    3-D  plot(vector  drawing  method). 

Consequently,  with  the  consideration  of  all  the  constraints  mentioned  above, 
another  type  of  format  was  made  to  see  the  details  of  the  beam  pattern.  Now,  a  plot 
was  made  using  the  function  SURMAT  about  the  6-^  plane  instead  of  the  x-y  plane. 
This  eliminated  the  previous  multiple  value  problem.  The  results  showed  the  most 
desirable  presentation  among  all  those  tried.  As  shown  in  Figure  3.7  it  gives  a  clearer 
appreciation  of  the  directivity  pattern  in  conjunction  with  0  and  <I>.  The  program  is 
hsted  in  Appendix  E. 

We  can  relate  this  three-dimensional  graph  (Figure  3.7)  to  the  two-dimensional 
graphs  (Figures  3.3  and  3.4)  for  interpretation.    As  shown  in  Figure  3.3  the  x-axis  is 
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the  angle  and  the  y-axis  is  the  directivity  which  is  expressed  in  dB.  Figure  3.4  is  the 
circular  graph  with  angle  working  counterclockwise  in  10^  increments  and  each  circle  is 
an  increment  of  10  dB  with  the  center  of  -  40  dB.  Figure  3.7  is  the  three-dimensional 
graph  from  the  view  point  of  20°  from  the  x-z  plane,  and  60°  from  the  x-y  plane.  Here, 
the  X-axis  is  O,  the  y-axis  is  0  and  the  z-axis  is  the  beam  pattern  (dB).  From  Figure 
3.3  four  major  points  can  be  found  at  approximately  50°,  130°,  230°  and  310°  of  6 
with  O  fixed  at  126°.  From  Figure  3.7  the  same  feature  can  be  determined  at  50°  of  6, 
126°  of  ^  as  seen  in  Figure  3.3.  Since  Figure  3.3  was  made  by  holding  <I>  constant  at 
126°,  one  can  look  along  the  <!>  =  126°  line  in  Figure  3.7  and  note  the  null  at  0  = 
90°,  the  maximum  at  0  =  50°,  and  the  relaiive  maximum  at  0  =  0°.  Note  the  heavy 
lines  in  Figure  3.7. 


DIRECriVITY  PATTERN 


K=10.5 

Dl=  0.5 

D2=0.5 

D3=0.5 

D4=0.5 

VIEW  ANGI.E=+20,+60 

P2=0  P3=0  P4=0 

A=L  B=l. 

C=l. 

Figure  3.7    3-D  plot(SURMAT  in  Q-<t>  plane). 
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From  the  circular  graph  (Figure  3.4)  four  major  points  can  be  found  approximatly  at 
40°,  140°,  220°  and  320°  of  <D  with  9  fixed  at  35°.  In  the  three-dimensional  graph  one 
of  these  major  points  is  seen  approximately  at  35°  of  G,  40°  of  <I>  and  the  other  major 
point  is  seen  approximately  at  35°  of  G,  140°  of  O.  Comparison  is  made  easier  when 
we  see  Figure  3.8  and  Figure  3.9  which  were  drawn  with  the  surfaces  cut  at  <P  =  126° 
and  G  =  35°  respectively.  Note  the  shift  in  axes  between  the  two  figures.  All  these 
points  of  two-dimensional  graphs  correspond  to  those  of  the  three-dimensional  graph. 


DIRECTIVITY  PATTERN 


K=10.5 

Dl=  0.5 

D2=0.5 

D3=0.5 

D4=0.5 

VIEW  ANGL£=+20,+60 

P2=0  P3=0   P4=0 

A=L  B=l. 

C=L 

Figure  3.8     3-D  plot(SURMAT  in  G-d)  plane). 

Inversely,  every  point  of  the  three-dimensional  graph  can  be  read  in  two-dimensional 
graphs  to  see  detailed  values  whenever  needed. 

Finally,  a  trial  was  made  to  draw  a  constant  value  plot  of  the  beam  pattern.  The 
first  approach  was  to  derive  the  functions  of  O  =  f  (G)  or  G  =  f  (O)  and  plot  these  on 
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Figure  3.9    3-D  plot(SURMAT  in  G-O  plane). 

the  G-O  plane.  However,  failing  to  derive  these  relationships  for  the  general  four 
element  case  made  this  method  of  making  a  contour  plot  of  a  fixed  F  value  impossible. 
Instead,  a  modified  version  of  the  existing  0-^  program  was  used.  It  sorted  on  a  small 
range  of  values  about  the  desired  value  of  F  to  find  the  corresponding  6-0  coordinates 
which  were  then  plotted  as  in  Figure  3.10.  This  type  of  plot  is  included  throughout  the 
study. 

With  these  forms  of  display  and  the  speed  of  the  FORTRAN  version  of  this 
program  we  can  now  use  this  method  to  study  the  directional  and  frequency  properties 
of  shaded  and  phased  simple  arrays. 
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TRACK  (-5  DB) 


Z  =  -5  DBK 

K=17.568 

D1=0  3  D2=0.3 

D3=0.15  D4=0. 

P2=  0   P3=180  P4=  0 

A=l.  B=l.  C=0. 


180.0 


PHI(DEG) 


Figure  3.10    contour  plot  (TRACK). 
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IV.  COMPARISON  WITH  SAMPLE  CASES 

In  making  comparisons,  some  simple  parameters  were  used  to  make  the 
comparisons  easier.  There  are  three  cases  included  in  this  chapter.  Each  case  has  a 
unique  set  of  input  data,  and  computer  generated  and  experimental  plots.  In  all  cases, 
a  two  points  were  picked  and  hand  calculations  made  to  verify  the  computer  output. 
An  example  of  the  hand  calculations  is  given  for  the  first  case  only.  Comparisons  were 
also  made  between  hand  or  computer  generated  values  and  the  two-dimensional  polar 
plot,  between  the  theoretical  two-dimensional  plot  and  the  experimental  data  plot, 
between  the  two-dimensional  plot  and  the  three-dimensional  plot,  and  fmally  between  a 
fixed  value  of  the  beam  pattern  (a  "track")  and  its  three-dimensional  plot. 

While  in  general,  9(0)  or  <P(6)  for  a  constant  F  could  not  be  obtained 
analytically,  in  a  few  cases,  included  below,  it  was  possible  to  obtain  the  functional 
dependencies.  These  are  discussed  in  the  following  cases. 

A.       CASE  1 

When  kD^  =  kD^  =  kD^  =  kD^  =  5.275, 

(P2  =  0,  (P3  =  (P4  =  :r, 

A  =  B  =  C  =  1 
equation  (2.37)  can  be  simplified  as 


H(e,(I))  =  lcos(5.275  sinGsinO)  -  cos(5.275  sin0cosO)|.  (4.1) 


Using  the  above  input  data  with  values  of  0  =  40°  and  O  =  0°  we  get  H(0,<I))  = 
0.984566  and  F  (0,O)  =  -0.135104  from  the  computer  output  as  seen  in  Appendix  G. 
If  we  use  the  above  simplified  equation  and  hand-computation,  we  get  H(0,C))  = 
0.98457  and  F(0,a))  =  -0.135068  . 

In  the  same  way,  if  9  =  90°  and  (P  =  0°  then  we  get  H  =  0.233302,  and  F  = 
-12.641622  from  the  computer  output  and  H  =  0.23330,  and  F  =  -12.6417  from  the 
hand-calculation.  These  hand-calculated  values  are  almost  the  same  the  the  computer 
output.  Comparisons  can  also  be  made  between  these  selected  values  and  the  two- 
dimensional    plot,    Figure    4.1a.      Additionally,    the    above    two-dimensional    plot 
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corresponds  well  in  shape  with  Figure  4.1b  which  is  from  LCDR  John  Butler's 
experiment  [Ref.  10]. 

If  we  compare  the  shape  of  the  two-dimensional  plot  (Figure  4.1a)  with  that  of 
the  three-dimensional  plot,  Figure  4.1c,  we  see  that  beam  pattern  of  the  two- 
dimensional  graph  is  identical  in  the  (l>  =  0^  plane  of  the  three-dimensional  plot. 


THETA  VS.  DIRECnVITY 


PHI=0  DEGREE 
K=10.55 
Dl=0.5  D2=0.5 
D3=0.5   D4=0.5 
P2=0  P3=180P4=180 
A=l.  B=l.  C=l. 


Figure  4.1a     2-D  plot  from  the  computer  output. 

A  two-dimensional  plot  of  equation  (4.1),  of  F  values  between  -10.9  dB  and  -10.0 
dB  results  in  Figure  4. Id.   For  F  values  less  than  -90.0  dB  we  get  Figure  4.1e. 
In  the  case  where  F  <  -90  dB,  H  <   10'^^''^^  which  is  small  enough  to  be  assumed  to 
be  0,  |cj  ~  Ic^l  or  e  =  nTi  n  =  0,1,2--,  and  equation  (2.40)  can  be  tan<I)  =   1  or  -1. 
From  these  resulting  equations  we  can  extract  <I>  =  45°  +  n7t/2,  n  =  1,2  — .   As  seen 
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in  Figure  4.1e,  when  F    <    -90  dB  the  asymptote  matches  the  nulls  in  the  0-<I) 
relationship  of  Figure  4.1c.  Note  the  nulls  at  9  =  0*^,  180°  for  all  <D,  and  <t>  =  45°. 
These  figures   of  constant  dB   level  aid  in  the  interpretation  of  Figure  4.1c.    The 
modified  P3D  program  (TRACE)  is  included  in  Appendix  F. 
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Figure  4.1b    2-D  plot  from  the  experiment. 
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DIRECnVITY  PATTERN 


K=10.55 

Dl=  0.5    D2=  0.5 
D3=  0.5     D4=  0.5 
VIEW  ANGLE=+20,+60 
P2=  0  P3=180  P4=180 
A=l.  B=l.  Ol. 


Figure  4.1c     3-D  plot  from  the  computer  output. 
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TRACK  (-10.  DB) 


Z  =  -10  DB 
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A=l.  B=l.  C=l. 
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Figure  4. Id    Track  (-10  dB). 
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TRACK  (-90  DB) 
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Figure  4.1e    Track  (-90  dB). 
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B.       CASE  2 

In  this  case,  a  triplet  array  is  investigated.  Referring  to  Figure  2.3,  the  fourth 
element  is  turned  off  and  its  opposing  element  is  moved  toward  the  center  of  the  array 
so  its  distances  is  half  that  of  the  two  elements.  Also,  the  source  2  is  180^  out  of  phase 
with  the  other  two  elements. 
When  kDj  =  kD^  =  5.2758, 

kDj  =  2.6379,  kD^  =  0, 

(P2  =  7t,  (P3  =  (p^  =  0, 

A  =  B=1,C  =  0 
the  equation  (2.37)  can  be  simplified  as 

H(0,<I>)  ^ {[cos(2.6379  sinGcosO)]^  +  [2  sin(5.2758  sinGsinO) 

-  sin(2.6379  sinGcosO)]^}  Vz.  (4.2) 

If  values  are  calculated  from  the  equation  (4.2)  as  in  case  1,  then  comparisons  can  be 
made  between  these  calculated  values  and  the  two-dimensional  plot  as  seen  in  Figure 
4.2a  (for  G  =  0^,  20°  and  90°,  and  <!>  =  90°,  we  get  F  =  -9.5424,  -2.7431  and  -3.6776 
respectively  from  hand-computation.).  Additionally,  the  two-dimensional  plot 
corresponds  well  in  shape  with  Figure  4.2b  from  the  experiment  (Ignore  the  ordering  of 
the  relative  phases  in  this  figures;  it  is  the  consequence  of  a  different  labeling 
convention). 

If  we  compare  the  shape  of  a  two-dimensional  plot  (Figure  4.2a)  with  that  of  the 
three-dimensional  plot  (Figure  4.2c),  we  see  that  beam  pattern  of  the  two-dimensional 
graph  is  identical  in  the  0  =  90°  plane  of  the  three-dimensional  plot. 

For  a  contour  plot  of  this  case,  a  two-dimensioal  plot  of  equation  (4.2)  of  F 
values  between  -5.9  dB  and  -5.0  dB  results  Figure  4. 2d.  For  F  values  between  -10.9  dB 
and  -10.0  dB,  we  get  Figure  4.2e.  For  F  values  less  than  -35.0  dB,  we  get  Figure  4.2f 
From  equation  (4.2),  we  get  F  =  -35.876  for  G  =  143°  and  (P  =  10°  from  hand- 
calculation.  This  point  is  seen  in  Figure  4.2f  Successive  examination  of  Figure  4. 2d 
through  4.2f  shows  how  the  contour  lines  migrate  as  lower  and  lower  values  of  F  are 
examined.  These  figures  of  constant  dB  level  aid  in  the  interpretation  of  Figure  4.2c. 
The  modified  P3D  program  (TRACE)  is  included  in  Appendix  F. 
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THETA  VS.  DIRECTrVITY 


PH1=  90  DEGREE 
K=17.586 
Dl=0.3  D2=0.3 
D3=0.15  D4=0. 
P2=180  P3=  0  P4=  0 
A=l.  B=l.  OO. 


Figure  4.2a     2-D  plot  from  the  computer  output. 
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TRIPOLE   BEAM   PATTERN    (dB) 
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Figure  4.2b    2-D  plot  from  the  experiment. 
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DIRECTIVITY  PATTERN 


K=17586 
Dl=  0.3     D2=  0.3 
D3=  0.15   D4=  0. 
VIEW  ANGLE=+20.+60 
P2=18GP3=0  P4=0 
A=l.  B=l.  C=0. 


Figure  4.2c    3-D  plot  from  the  computer  output. 
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TRACK  (-5.0  DB) 


Z  =  -5.0  DB 
K=17.586 
Dl=0.3   D2=0.3 
D3=0.15  D4=0. 
P2=160P3=0   P4=0 
A=l.  B=l.  C=0. 


60.0  100.0 

THETA(DEG) 
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140.0 
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180.0 


Figure  4.2d    Track  (-5  dB). 
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TRACK  (-10.0  DB) 


Z  =  -10.0  DB 

K=17.586 
Dl=0.3   D2=0.3 
D3=0.15  D4=0. 
P2=180  P3=  0   P4=  0 
A=L  B=L  C=0. 


40.0  aO.O  80.0  100.0  120.0  140.0  160.0  180.0 


THKrA(DEG) 


Figure  4.2e    Track  (-10  dB). 
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'i'KACK  (-35.0  DB) 

Z  =  -35.0  DB 
K=17.5a6 
Dl=0.3   D2=0.3 
D3=0.15  D4=0. 
P2=180P3=0   P4=0 
A=l.  B=l.  C=0. 

PHI(DEG) 

0          10  0          80.0          30.0          40  0          60.0          60.0          70  0          80  0          9C 

0  0             20  0             40.0             60.0             80.0            100.0            120.0            140.0            160.0            Ifi 

THhTA(DEG) 

0.0 

Figure  4.2f    Track  (-35  dB). 
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C.       CASE  3 

In  this  case,  a  doublet  array  is  investigated.  Referring  to  Figure  4.2a,  this  is  the 
same  as  the  first  case  except  that  the  third  and  fourth  elements  are  turned  off  and  all 
elements  are  in  phase. 
For      kD^  =  kD.  =    kDj  =  kD^  =  5.275, 

A=1,B  =  C  =  0 

equation  (2.37)  can  be  simplified  as 

H(e,<I))  =  lcos(5.275  sinesin<D)| 

As  in  cases  1  and  2,  for  9  =  0°,  35°,  90°  and  O  =  90°,  we  get  F  =  0,  -0.058  and 
-5.459  respectively  from  hand-computation.  These  points  are  found  in  Figure  4.3a  and, 
the  two-dimensional  plot  again  corresponds  well  in  shape  with  Figure  4,3b  from  the 
experiment. 

If  we  compare  the  shape  of  a  two-dimensional  plot  (Figure  4.3a)  with  that  of  the 
three-dimensional  plot  (Figure  4.3c  below),  we  see  that  the  beam  pattern  of  the  two- 
dimensional  graph  is  identical  in  the  ^  =  0°  plane  of  the  three-dimensional  plot. 

For  a  constant  value  of  H(0,<I>),  we  can  follow  the  same  procedure  to  get  the 
equation  for  contour  plots  as  in  case  1  and  2.  The  resultant  figures  are  Figures  4.3d, 
4.3e  and  4.3f  Figure  4.3d  is  drawn  for  the  F  values  between  -5.9  dB  and  -5.0  dB  and 
Figure  4.3e  is  for  the  F  values  between  -8.9  dB  and  -8.0  dB.  Figure  4.3f  is  for  F  values 
less  than  -35.0  dB. 
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THETA  VS.  DIRECTIVITY 


PH]=  90  DEGREE 
K=10.55 
Dl=0.5  D2=0.5 
D3=0.5   D4=0.5 
P2=0  P3=0  P4=0 
A=L  B=0.  C=0. 


Figure  4.3a    2-D  plot  from  the  computer  output. 
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DIPOLE   BEAM   PATTERN    (dB) 
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Figure  4.3b    2-D  plot  from  the  experiment. 
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DIRECriVITY  PATTERN 


K=10.55 

Dl=0.5    D2=0.5 
D3=  05    D4=  0.5 
VIEW'  ANGLE=+20,+60 
P2=0  P3=0   P4=0 
A=L  B=0.  C=0. 


Figure  4.3c    3-D  plot  from  the  computer  output. 
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TRACK  (-  5.0  DB) 


Z  =  -5.0  DB 

K=10.55 
Dl=0.5  D2=0.5 
D3=0.5   D4=0.5 
P2=0   P3=0  P4=0 
A=L  B=0.  0=0. 
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Figure  4.3d    Track  (-5  dB). 
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TRACK  (-8  DB) 


0.0  ZOS)  40JD  60.0  80.0  100.0  120.0 

THETA(DEG) 


Z  =  -8.0  DB 
K=10.55 
Dl=0.5  D2=0.5 
D3=0.5   D4=0.5 
P2=  0   P3=  0  P4=  0 
A=l.  B=0.  C=0. 


140.0  160.0  180.0 


Figure  4.3e     Track  (-8  dB). 
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TRACK  (-35.0  DB) 
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Z  =  -35.0  DB 
K=10.55 
Dl=0.5   D2=0.5 
D3=0.5   D4=0.5 
P2=  0   P3=  0  P4=  0 
A=L  B=0.  C=0. 


0.0  20.0  40.0  80.0  80.0  lOO.O  120.0  140.0  160.0  180.0 

THETA(DEG) 


Figure  4.3f    Track  (-35  dB). 
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V.  EVALUATION 

In  the  last  chapter,  we  looked  at  three  simple  cases.  These  established  the 
accuracy  of  the  five  computer  programs.  For  example,  in  comparisons  between  hand 
calculation  and  computer  output,  both  sets  of  values  were  proven  nearly  identical. 
Again,  due  to  the  complexity  of  calculations,  not  all  values  were  hand-calculated. 

When  a  two-dimensional  circular  graph  was  compared  with  the  selected  hand- 
computed  values,  both  groups  of  values  were  again  virtually  identical.  As  another  test 
of  the  validity  of  the  five  computer  programs,  some  of  the  computer  generated  two- 
dimensional  plots  were  compared  against  those  from  an  experiment  [Ref  9:  p.l8|. 
Visual  comparisons  of  the  two  types  of  plots  were  similar  in  shape,  but  values  in  dB 
were  somewhat  diiferent.  These  differences  might  be  caused  by  the  normalization  of 
the  angular  dependance  of  the  sound  pressure.  Also,  because  of  the  directivity  of  the 
real  sources,  there  is  a  slight  decrease  in  dB  level  of  the  experimental  plots  when  the 
array  pointed  away  from  the  receiver.  However,  they  are  small  enough  that  the 
comparisons  still  seem  quite  reasonable. 

Next,  two-dimensional  versus  three-dimensional  graphs  were  compared. 
Agreement  here  shows  the  programs  are  working  properly.  By  varying  the  range  of  8 
and  <I)  that  are  plotted,  and  also  by  varing  the  view  angles,  values  and  the  shape  of  the 
curve  can  be  visualized  without  great  difficulty. 

Tracks  of  certain  directivity  values  were  drawn  using  the  modified  P3D  program 
(TRACE)  and  were  compared  to  the  three-dimensional  plots.  The  shape  of  the 
resulting  tracks  corresponded  well  with  those  of  the  three-dimensional  patterns,  again 
showing  proper  functioning  of  the  programs. 

With  all  the  comparisons  made  in  Chapter  IV,  the  five  computer  programs 
developed  in  this  thesis  were  proven  effective  in  predicting  the  directivity  pattern  in  all 
directions  of  space  with  or  without  differences  in  the  strength  of  the  sources,  inter- 
source  distances  and  phases. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  following  conclusions  are  offered. 

•  The  P2DICIR,  P2DPCIR,  P2DIFIX,  P2DPFIX  programs  successfully 
predict  the  directivity  beam  pattern  in  two  dimensions. 

•  The  P3D  program  predicts  and  visualizes  the  directivity  beam  pattern  in  the 
three-dimension  space. 

•  They  allow  the  presentation  of  complete  beam  patterns,  which  contain 
considerable  information,  in  a  succinct,  easy  to  comprehend  and  interpret  form. 

B.  RECOMMENDATIONS 

The  following  recommendations  and  areas  of  further  investigation  is  suggested. 

•  Further  studies  for  computer  programs  to  visualize  the  beam  pattern  in  the 
x-y-z  space  are  needed  to  eliminate  hidden  lines  in  three-dimensional  plotting. 

•  Transcendental  equations  need  to  be  solved  to  get  the  contour  plot  for  fL\ed 
F  values. 
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APPENDIX  A 
BASIC  PROGRAMS  (HP  COMPUTER) 


IIKMKXMXMMWXMXKXXXIOCKXKM        INFORMATION        XXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXM 

*  THE  OBJECTIVE  OF  THESE  PROGRAMS  IS  TO  CALCULATE  AND  » 

»        DRAW  THE  DIRECTIVITY  PATTERN  IN  TWO-DIMENSION.  » 

»        THIS  IS  WRITTEN  IN  BASIC  LANGUAGE  AND  PLOTTER  IS  HP7090A.      » 
»  » 

MXXMMKXKXKXxxxxxxxxxxw  VARIABLE  DEFINITION  xxxxxxxxkxxkkxxmxxxmxkkxxxx 

»  Z  :  DIRECTIVITY  IN  dB  » 

»  P2  :  PHASE  (DEGREE)  OF  SOURCE  A  » 

»  P3  :  PHASE  (DEGREE)  OF  SOURCE  B  * 

»  P<¥    :  PHASE  (DEGREE)  OF  SOURCE  C  » 

«  P5  :  VARIABLE  TO  CHANGE  DEGREE  TO  RADIAN  « 

»  P6  :  RADIAL  VALUE  OF  PI  » 

»  P7  :  RADIAL  VALUE  OF  P2  » 

»  PS  :  RADIAL  VALUE  OF  P5  » 

«  K  :  WAVE  NUMBER  « 

»  D  :  INTER-SOURCE  DISTANCE  » 

»  P  :  PHI  IN  DEGREE  » 

»  I  :  THETA  IN  DEGREE  » 

»  PI  :  3.1^15917  * 

»  « 

MKMMKXXXMXXXMKXXXXXKXXKMXXXXMXMMXXXXXMKKMKXXXKXXXMXXXMMXXXXXKKXXXXXXKXKM 
XXXXKXXXXXXKXXKXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

HP  PROGRAM  FOR  THE  RECTANGULAR  GRAPH 

XXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXX 


INPUT  PARAMETERS 

5  OPTION  BASE  1 

6  RAO 

10  DIM  Z(360) 

20  DISP  "WHAT  IS  THE  VALUE  OF  KD?  KD=   "  !  KD=K»D 

30  INPUT  KD 

60  DISP  "WHAT  IS  THE  VALUE  OF  P2(DEG)?  P2=  "!        P2=PHASE  OF  A 

61  INPUT  P2 

62  DISP  "WHAT  IS  THE  VALUE  OF  P3(DEG)?  P3=  "  »       P3=PHASE  OF  B 

63  INPUT  P3 

6^  DISP  "WHAT  IS  THE  VALUE  OF  P3(DEG)?  P<*=  "  !       P^=PHASE  OF  C 

65  INPUT  P4 

66  DISP  "WHAT  IS  THE  VALUE  OF  A?  A=   " 

67  INPUT  A 

68  DISP  "WHAT  IS  THE  VALUE  OF  B?  B=   " 
6  9  INPUT  B 

70  DISP  "WHAT  IS  THE  VALUE  OF  C?  C=   " 

71  INPUT  C 

72  DISP  "IS  P  FIXED?  (YES>AN=1  ,  NOjAN=0)   AN=   " 

73  INPUT  AN 

7<t  IF  AN=1  THEN  DISP  "WHAT  IS  THE  VALUE  OF  P(DEG)?  P=  "ELSE  GOTO  76 

75  INPUT  P 

76  DISP  "IS  I  FIXED?  (YES>ANS=1  ,  NOjANS=0)   ANS=  " 

77  INPUT  ANS 

78  IF  ANS=1  THEN  DISP  "WHAT  IS  THE  VALUE  OF  KDEG)?  1=  "ELSE  GOTO  82 

79  INPUT  I 

83  IF  AN=0  AND  ANS=0  THEN  GOTO  86 

84  IF  AN=0  AND  ANS=1  THEN  GOTO  90 

85  IF  AN=1  AND  ANS=0  THEN  GOTO  93 

86  DISP  "»">"K*0=">KD,"  ","»" 

87  DISP  "*PHASE  OF  A( DEG )=" )P2>"PHASE  OF  B( DEG )=" jP3>"PHASE  OF  C( DEG)=")P<i,"»" 

88  DISP  "»"v"A=">A,"  B="lB,"C="lC,"»"  3  GOTO  99 

90  DISP  "*")"K*D="-,KD,"   I(DEG)=")I,"  ","*•' 

91  DISP  "JfPHASE  OF  A( DEG  )=" >P2j"PHASE  OF  B( DEG  )="jP3>"PHASE  OF  C( DEG )="jP4,"»" 

92  DISP  "*")"A=";A,"  B=">B,"C=" )C,"*"  3  GOTO  103 

93  DISP  "»"j"KJ*D="vKD,"   I(DEG)=">I,"  ","*" 
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9<*  DISP  "»PHASE  OF  A( DEG  )=">P2)"PHASE  OF  B( DEG  )="jP3j"PHASE  OF  C( 0EG)=")P4,"»" 
95  OISP  "*"i"A=";A,"  8="iB,"C=" jC,"*"  3  GOTO  107 

99  OISP  "XMKKXXXMXKXKKMKXKKXXXKXXXXKKXKMXKXXKXXXMXXXXMXXKXMMXX" 

100  DISP  "  " 

101  DISP  "   "}"  THETA"J"    PHI","  DIRECTIVITY" 

102  FOR  P=l  TO  360  3  GOTO  110 

103  DISP  "XXXXXKXXXXXXXXXKXXXXXXKXXXXXXXXWXXXXXXXXXXXXKXXKKXXKK" 

10<f  DISP  "  " 

105  DISP  "   "J"  THETA">"     PHI","  DIRECTIVITY" 

106  FOR  P=l  TO  360  3  GOTO  IIS 

107  DISP  "XXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXKKXKXXKXXXXXKX" 

108  DISP  "  " 

109  DISP  "  ";•'   THETA">"     PHI","  DIRECTIVITY" 

MAIN  PROGRAM 

110  FOR  1=1  TO  360 

115  P5=PI  /ISO 

116  P6=P2»PS  3  P7=P3*PS  3  P8=P4»P5 

117  Q=P-90 

118  Q1=SIN  (Q*PS) 

119  I1=SIN  (I*P5) 

120  Q2=SIN  (P*PS) 

121  Wl=KD*Il*Q2/2 
121  H2=KD*Il*Ql/2 

127  W=COS  (W1)+A*C0S  ( W1+P6 )+B»COS  ( W2-P7 )+C»COS  (W2+P8) 

128  G=SIN  (W1)-A*SIN  ( W1+P6 )+B»SIN  (M2-P7)-C»SIN  (W2+P8) 
130  H=SQR  (W**2+GJ«»2)/(1+A+B+C) 

140  Z(P)=20*LGT  (H) 

143  T=20*LGT  (ABS  (SIN  (.S*KD»I1))) 

ISO  DISP  USING  1S3  >  I,P,Z(P) 

153  IMAGE  5X,3D.D,5X,3D.D,1SX,SDDD.DDD,15X,SD0D.DD0 

155  !   IF  Z(I]<-20  THEN  I=I+.l  3  GOTO  119 

156  !   1=  IP(I) 

157  !   IF  Z(P)<-20  THEN  P=P+.l  3  GOTO  117 

158  !   P=  IP(P)) 

160  IF  AN=1  AND  ANS=0  THEN  GOTO  163  ELSE  GOTO  161 

161  IF  AN=0  AND  ANS=1  THEN  GOTO  164  ELSE  GOTO  163 

162  NEXT  P  3  GOTO  96 

163  !   NEXT  I 

164  NEXT  P 


PLOTTING  WITH  HP7090A 


170 

PLOTTER  IS  70S 

180 

CLEAR 

190 

GCLEAR 

200 

GRAME 

210 

LOCATE  15,145,10,80 

220 

SCALE  0,390,-240,0 

230 

CLIP  0,360,0,-240 

240 

AXES  5,5,0,-240,10,10,2 

250 

LINE  TYPE  1 

260 

GRID  20,-20 

270 

FRAME 

280 

LORG  6 

290 

CSIZE  4 

300 

FOR  1=0  TO  360  STEP  40 

310 

MOVE  1,-241 

320 

LABEL  I 

330 

NEXT  I 

340 

LORG  8 

350 

FOR  I =-240  TO  0  STEP  40 

360 

MOVE  -1,1 

370 

LABEL  I 

380 

NEXT  I 

390 

CSIZE  5 

400 

LORG  5 

410 

MOVE  180,40 

420 

CSIZE  * 

430 

PEN  2 

440 

LABEL  "PHI  VS  DIRECTIVITY" 

450 

LORG  5 

460 

MOVE  180,20 

470 

CSIZE  4 

480 

LABEL  "(THETA  =  45  DEGREE)" 

490 

CSIZE  5 

500 

LORG  5 

510 

MOVE  180,-26- 

520 

CSIZE  4 

52 


530  LABEL  "PHKDEG)" 
540  LDIR  90 
550  MOVE  -20,-120 
560  DEG 
570  LDIR  90 
580  CSIZE  4 

590  LABEL  "DIRECTIVITY(dB )" 
600  CSIZE  5 
610  MOVE  0,0 
620  FOR  P=l  TO  360 
630     DRAW  P,Z(P) 
640  NEXT  P 
650  PEN  0 
700  END 
117392 
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HP  PROGRAM  FOR  CIRCULAR  GRAPH 


INPUT  PARAMETERS 


KD=K»D 

P2=PHASE  OF  A 
P3=PHASE  OF  B 
P4=PHASE  OF  C 


5  OPTION  BASE  1 

6  RAD 

10  DIM  X(360) 

11  DIM  Y(360) 

12  DIM  Z(360) 

20  DISP  "WHAT  IS  THE  VALUE  OF  KD?   KD=   "  » 
30  INPUT  KD 

60  DISP  "WHAT  IS  THE  VALUE  OF  P2(DEG)?  P2=  "! 

61  INPUT  P2 

62  DISP  "WHAT  IS  THE  VALUE  OF  P3(DEG)?  P3=  "  ! 

63  INPUT  P3 

64  DISP  "WHAT  IS  THE  VALUE  OF  P4rDE6)?  P4=  "  ! 

65  INPUT  P4 

66  DISP  "WHAT  IS  THE  VALUE  OF  A?   A=   " 

67  INPUT  A 

68  DISP  "WHAT  IS  THE  VALUE  OF  B?   B=   " 

69  INPUT  B 

70  DISP  "WHAT  IS  THE  VALUE  OF  C?   C=   " 

71  INPUT  C 

72  DISP  "IS  P  FIXED?  (YESvAN=l  ,  NOjAN=0)   AN=   " 

73  INPUT  AN 

74  IF  AN=1  THEN  DISP  "WHAT  IS  THE  VALUE  OF  P(DEG)?  P=  "ELSE  GOTO  76 

75  INPUT  P 

76  DISP  "IS  I  FIXED?  (YESjANS=l  ,  NOjANS=0)   ANS=  " 

77  INPUT  ANS 

78  IF  ANS=1  THEN  DISP  "WHAT  IS  THE  VALUE  OF  I(DEG)?  1=  "ELSE  GOTO  82 

79  INPUT  I 

82  DISP  "XKKXXMMKXXXKXMKXXKXXXXXXKXXXXMXXMXXMXXXKKXXXKXKXKXXXK" 

83  IF  AN=0  AND  ANS=0  THEN  GOTO  86 

84  IF  AN=0  AND  ANS=1  THEN  GOTO  90 

85  IF  AN=1  AND  ANS=0  THEN  GOTO  93 

86  DISP  "*"r*K*D="jKD,"  "  ,"*" 

87  DISP  "«PHASE  OF  A( DEG  )=")P2 j"PHASE  OF  B( DEG  )=" jP3)"PHASE  OF  C( DEG )=">P4,"»" 

88  DISP  "*">"A="vA,"  B="iB,"C="jC,"»"  3  GOTO  99 

90  DISP  "*"j"K*D="jKD,"   I(DEG)=">I,"  ","*" 

91  DISP  "*PHASE  OF  A( DEG  )  =  " iP2j"PHASE  OF  B( DEG  )=" }P3r'PHASE  OF  C( DEG )=" jP4,"*" 

92  DISP  "*")"A=")A,"  B="jB,"C="}C,"«"  3  GOTO  103 

93  DISP  "*"J"K«D="JKD,"   I(DEG)="jI,"  ","»" 

94  DISP  "*PHASE  OF  A( DEG  )=" iP2 j"PHASE  OF  B( DEG)=")P3 J"PHASE  OF  C( DEG)=">P4,"»" 

95  DISP  "»"}"A="jA,"  B=">B,"C=">C,"»"  S  GOTO  107 

99  DISP  "KXKXXXXXXXXMXKXXXXXXXXKXXXXKXXXXXXXXKXXXKKXXKXXKKKXXK" 

100  DISP  "  " 

101  DISP  "   "5"  THETA";"    PHI","  DIRECTIVITY" 

102  FOR  P=l  TO  360  3  GOTO  110 

103  DISP  "KXXXXKKXXKXXKXKKXKXXKKXKKKXKXXXXKXKKXXKKXKXKKXKXKKKMM" 

104  DISP  "  " 


105  DISP 


i"   THETA"j" 


PHI' 


DIRECTIVITY' 


106  FOR  P=l  TO  360  3  GOTO  115 

107  DISP  "XXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX" 

108  DISP  "  " 

109  DISP  "   "}"  THETA"j"    PHI","  DIRECTIVITY" 


MAIN  PROGRAM 


110  FOR  1=1  TO  360 
115  P5=PI  /180 
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116  P6=P2»P5  3  P7=P3«PS  S  P8=P^«PS 

117  Q=P-90 

118  Q1=SIN  (Q«PS) 

119  I1=SIN  (I*PS) 

120  Q2=SIN  (P«P5) 

121  Wl=KD*Il*Q2/2 
121  W2=KD*Il*Ql/2 

127  W=COS  (W1)+A*C0S  ( W1+P6  )+B*COS  ( W2-P7)+C»C0S  (W2+P8) 

128  G=SIN  (W1)-A*SIN  ( W1+P6 )+B*SIN  ( W2-P7)-C»SIN  (W2+P8) 
130  H=SQR  (W**2+G»*2)/(l+A+B+C) 

1^0   Z(I)  =  20*LGT    (H)  +  2'»0 

l<tl  X(I)=Z(I)*COS    (I»P5) 

l<+2   Y(IJ=Z(I)*SIN    (I*P5) 

1«   T=20»LGT    (ABS    (SIN    (.S»KD*I1))) 

ISO  DISP  USING  153  >  I,P,Z(P) 

153  IMAGE  5X,SDDD.00D,15X>SDDD.DD0 

155  !   IF  Z(I)<-20  THEN  1=1*. 1  3  GOTO  119 

156  !   1=  IP(I) 

157  !   IF  Z(P)<-20  THEN  P=P+.l  3  GOTO  117 

158  !   P=  IP(Pn 

160  IF  AN=1  AND  ANS=0  THEN  GOTO  163  ELSE  GOTO  161 

161  IF  AN=0  AND  ANS=1  THEN  GOTO  16^  ELSE  GOTO  163 

162  NEXT  P  3  GOTO  96 

163  NEXT  I 

164  !   NEXT  P 


PLOTTING  WITH  HP7090A 


170  PLOTTER  IS  70S 

180  CLEAR 

190  GCLEAR 

200  GRAME 

210  LOCATE  25,120,5,80 

220  SCALE  -260,320,-220,2«i0 

230  CLIP  -2'+0,2<+0,-2<t0,2<f0 

235  XAXIS  0,20,-240,240 

240  YAXIS  0,20,-240,240 

250  LINE  TYPE  1 

270  FRAME 

280  LORG  6 

290  CSIZE  4 

300  FOR  1=0  TO  360  STEP  40 

310     MOVE  1,-5 

320  LABEL  1-240 

321  MOVE  -1,1 

322  LABEL  1-240 
330  NEXT  I 

340  LORG  6 

341  MOVE  0,-S 

342  LABEL  -240 

350  FOR  I=-240  TO  -60  STEP  60 
360     MOVE  -1,1 

370  LABEL  -240-1 

371  MOVE  1,-5 

372  LABEL  -240-1 
375  NEXT  I 

380  MOVE  0,0 

381  FOR  1=0  TO  240  STEP  60 

382  FOR  J=0  TO  360 

383  X=I»COS  (J*PI  /ISO) 

384  Y=I*SIN  (J*PI  /180) 

385  DRAW  X,Y 

386  NEXT  J 

387  NEXT  I 
390  CSIZE  5 
400  LORG  5 
410  MOVE  0,310 
420  CSIZE  8 
430  PEN  2 

440  LABEL  "THETA  VS  DIRECTIVITY" 

450  LORG  SO 

460  MOVE  0,285 

470  CSIZE  4 

480  LABEL  "(PHI  =  90  DEGREE)" 

490  CSIZE  5 

500  LORG  5 

510  MOVE  300,0 

520  CSIZE  4 

530  LABEL  "O(DEG)" 

550  MOVE  -300,0 

580  CSIZE  4 

590  LABEL  "180(DEG)" 
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610  MOVE  0,0 
620  FOR  I-l  TO  360 
630     DRAW  X(IJ,Y(IJ 
6<t0  NEXT  I 
650  PEN  0 
700  END 
11117^ 
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APPENDIX  B 
P2DIFIX,  P2DPFIX,  P2DICIR  AND  P2DPCIR 

KKKXKKXXKKXKXXXXKKKXKKK   INFORMATION   XKXXKXXXXXXXXKXXXXXXXXXXKKKXXXKXKK 

*  M 

«  THE  OBJECTIVE  OF  THESE  PROGRAMS  IS  TO  CALCULATE  AND  » 

«  DRAW  THE  DIRECTIVITY  PATTERN  IN  TWO-DIMENSION.  « 

»  THE.  DRAWING  METHOD  IS  THE  SUBROLTTINE  CURVE  OF  DISSPLA.  » 

»  1(9)  AND  P(<P)  CAN  BE  USED  INTERCHANGEABLY  FOR  « 

»  P2DPFIX  AND  P2DICIR  PROGRAMS.  » 

»  » 
xxxxxKKKxxxxxxxxxKxxxx   VARIABLE  DEFINITION  xxxxxkxxxxxxxxkkxkxxkxxkxxx 

*  2  :  DIRECTIVITY  IN  dB  « 
«  PI  :  PHASE  (DEGREE)  OF  SOURCE  A  « 
»  P2  :  PHASE  (DEGREE)  OF  SOURCE  B  » 
«  P3  :  PHASE  (DEGREE)  OF  SOURCE  C  » 
»  P4  :  VARIABLE  TO  CHANGE  DEGREE  TO  RADIAN  » 
«  PS  :  RADIAL  VALUE  OF  PI  « 
«  P6  :  RADIAL  VALUE  OF  P2  » 
»  P7  :  RADIAL  VALUE  OF  P3  « 
»  K  :  WAVE  NUMBER  * 
»  Dl  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  1  « 
»  02  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  A  » 
«  03  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  B  « 
»  04  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  C  * 
»  P  :  PHI  IN  DEGREE  » 
»  I  :  THETA  IN  DEGREE  » 
«  PH  :  PHI  IN  RADIAN  » 
»  TH  :  THETA  IN  RADIAN  » 
«  PHI  :  DUMMY  VARIABLE  OF  PHI  » 
»  THETA  :  DUMMY  VARIABLE  OF  THETA  » 
«  PI  :  3.1415917  » 
»  » 

KXKMMWKKKXKKMKKKKXMKXK   VARIABLE  DECLARATION   KWXXXKKKKXXKXWKKXXKXKKKXWM 

C 

c 

KXKKXXXXKXXXXXXKXKXKXXXXXKXKXXKXXXXXXXXXXXMXXXXWXKXXKXXXKXXKKXXXXXXXXXKK 

PROGRAM  FOR  THE  TWO-DIMENSIONAL  RECTANGULAR  GRAPH 

C 

DIMENSION  PHI ( 361 )  ,THETA( 361 )  ,Z( 360  ) 

REAL  P4,P5>P6,P7,Wl,W2,W3>W4>W,G,H,K,01,02,03,D<f,PH,TH>PI,a,b,c 

INTEGER  P,P1,P2,P3,I 
C 
c 

C  INPUT  PARAMETERS 

C 

c 

PI=3.1'*1S917 
K=S.S 

01  =  1. 

02  =  1. 
03=1. 
04=1. 
Pl=  180 
P2=  0 
P3=  180 
A  =  1. 
B  =  1. 
C  =  1. 

I  =  25:  USE  P  FOR  P2DPFIX 
C 
c 

C  MAIN  PROGRAM 

C 

C 


P4=PI/180 
P5=P1*P4 
P6=P2*P4 
P7=P3*P4 

00  100  1=1,360:  MOVE  C  TO  LINE  NO.  76  FOR  P20PFIX 
00  200  P=l,360 
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TH=I*P<i 

PH=P*P4 

W1=K«D1«SIN( PH  )«SIN{ TH ) 

W2=K»D2*SIN( PH )«SIN( TH )_P5 

H3=K*D3*SIN( TH )*COS( PH )-P6 

W4=K*D4*SIN( TH )«COS( PH )+P7 

W=COS( Wl )+A*COS( W2 )+B»COS( HZ )+C»COS( H4 ) 

G=SIN(  Wl  )-A*SIN(  W2  )-B»(SIN(  W3  )+C*SIN(  W4 ) 

H=SQRT( W*W+G«G )/( 1+A+B+C ) 

2(P)=20»ALOG10(H) 

IF(Z(PJ.LE.-'+0.  )  THEN 
Z(P)=-40 

ELSE 

2tP)=Z(P) 

END  IF 

THETA(P)=I 

PHI(P)=P 

WRITE! 8,1)  PHI(P),2(P) 
1           FORMAT! 2X,2( ElO. 3, SX)) 
200      CONTINUE 
100   CONTINUE 
C 
C 

C         A  PROGRAM  FOR  PLOTTING  BY  SHERPAOR  BY  TEK618 

C 

C 

CALL  MEDBUF 

CALL  COMPRS 
C     CALL  TEK618:  MOVE  C  TO  DRAW  USING  TEK618 

CALL  PHYSORCl. ,1.  ) 

CALL  PAGE! IS. 5, 12.  ) 

CALL  AREA2D( 10.5,9.1) 

CALL  FRAME 

CALL  COMPLX 

CALL  HEIGHT! .2) 

CALL  XNAMEt  'PHI (DEG)   $',10) 

CALL  YNAME( 'DIRECTIVITY (DB)$', IS) 

CALL  YTICKS(8) 

CALL  XTICKS(8) 

CALL  GRAF(0,<fO,360,-40,<»,0) 

CALL  GRID! 2,2) 

CALL  CURVE! PHI ,Z, 359,0) 

CALL  MESSAG! •THETA=25  DEGREE ' ,15,11. ,8.5) 

CALL  MESSAG! 'KsS. 5' ,5,11. ,8.2) 

CALL  MESSAG! 'Dl=l.  02=1. ' ,11,11. ,7. 9) 

CALL  MESSAG! 'D3=l.  D^=l. ' ,11 ,11. ,7.6 ) 

CALL  MESSAG!  'P2=180  P3=0  P'+slBO  '  ,18,11.  ,7.3  ) 

CALL  MESSAG! 'A=l.  B=l.  C=l. ' ,14,11. ,7.0 ) 

CALL  RESET! 'HEIGHT' ) 

CALL  HEIGHT!  .<♦) 

CALL  MESSAG! 'PHI  VS.  DIRECTIVITY' ,19,3.1,9.8) 

CALL  RESET!  'ALL'  ) 

CALL  NOCHEK 

CALL  ENOGRIO) 

CALL  PHYSOR! 11.85,7.6) 

CALL  AREA2D!3.5,2.5) 

CALL  FRAME 

CALL  ENDPLiO) 

STOP 

END 
C 
C 

KKKKKKKMKKKXKKKMKKKMMKWKWKKKICKWICKMKXKMMKKKKKMKXKKMKKKKKKMICKKICKKKMMXMKXKM 

PROGRAM  FOR  THE  CIRCULAR  TWO-DIMENSIONAL  GRAPH 

KKKKKKKKXXKKKKKKK)CKKKXKKKMXKKXKK)CKKKMKK*KKKXKKKMKMKKMKMKKXMKKXMXKKKKKKKM 
C 

DIMENSION  PHI!  360  )  ,R!  360  )  ,THETAI  360  )  ,ZI  360  )  ,XI  4 )  ,YI  <* ) 

REAL  P4,P5,P6,P7,Wl,W2,W3,W<+,W,G,H,K,Dl,D2,03,D<f,PH,TH,PI,a,b,c 

INTEGER  P,P1,P2,P3,I,Q 

C 

c 

C  INPUT  PARAMETERS 

C 

c 

PI=3. 1415917 

K=5.5 

Dl=l. 

D2=l. 

D3=l. 

D4=l. 

Pl=  180 

P2=  0 

P3=  180 

A  =  1. 
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B  =  1. 
C  =  1. 
I    =   25:    USE   P   FOR  P20PCIR 


1 

200 

100 

c 
c— 

c 

c 

c 


MAIN  PROGRAM 


P<f=PI/180 
P5=Pl*P<t 

P7=P3*P<+ 

DO  100  1=0,260:  MOVE  C  TO  LINE  NO.  174  FOR  P2DICIR 
DO  200  P=0,360 

TH=I«P<t 

PH=P*P4 

H1=K*D1*SIN( PH )«COS(  TH ) 

W2=K*D2*SIN( PH )*COS(  TH )+P5 

W3=K*D3*SIN( TH )*COS( PH )-P6 

W«^=K»D4*SIN(  TH  )*COS(  PH  )+P7 

H=COS( HI J+A*COS( H2  )+B*COS( W3  )+C»COS( W4 ) 

G=SIN( Wl  )-A*SIN( W2 )+B*SIN{ W3 )-C*SIN( W4 ) 

H=SQRT( W*W+G*G )/( 1+A+B+C ) 

21P)=20*AL0G10(H) 

IF(2(P).LE.-<»0.  J  THEN 


2(P)=0. 
ELSE 

Z(P)='+0.+Z(P) 
END  IF 
THETA(P)=I 
PHI(P)=P 

WRITE(8,1)  PHI(P),Z(P) 
F0RMAT(2X,2(E10.3,SX)) 
CONTINUE 
CONTINUE 


A  PROGRAM  FOR  PLOTTING  BY  SHERPAOR  BY  TEK618 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
STOP 
END 


MEDBUF 

COMPRS 

TEK618:  MOVE  C  TO  THE  UPPER  LINE  HHEN  DRAWING  WITH  TEK618 

PHYSOR( 1.5,1.  ) 

PAGEf 15.5,12.  ) 

AREA2D(3. ,8.  ) 

COMPLX 

HEIGHT( .2) 

POLAR(  3.  l<t  15917/180., 10., 'i.,<i.  ) 

GRID(1,1) 

CURVE( PHI ,Z, 360,0) 

MESSAG( 'THETA=25  DEGREE ' ,15,9. 7,7.7) 

MESSAGt  'K=5.5',5,9.7,7.4) 

MESSAG( 'Dl=l.  02=1. ',11,9.7,7.1) 

MESSAG(  'D3  =  l.  04=1. ' ,11,9. 7,6 .8 ) 

MESSAGt  •P2=180  P3=0  P4= 180 ' ,18,9. 7,6.5 ) 

MESSAGt  'A=l.  B=l.  C=l. ' ,14,9. 7,6 . 2 ) 

MESSAG(  '0'  ,1,3.9,8.1) 

MESSAG( '-10' ,3,3.7,7.1) 

MESSAGl  '-20'  ,3,3.7,6.1) 

MESSAG( '-30', 3, 3. 7, 5.1) 

MESSAG( '-40' ,3,3.7,4.1) 

MESSAGt  '0(DEG)',6,7.6,3.9) 

MESSAGt  '180' ,3,0., 3. 9) 

RESET(  'HEIGHT'  ) 

HEIGHT! .4) 

MESSAGt  'PHI  VS.  DIRECTIVITY* ,19,3.0,9.5) 

RESETt  'ALL'  ) 

NOCHEK 

ENDGR(O) 

PHYSORt 10. 95,6.8) 

AREA2D(3.5,2.5) 

FRAME 

DONE PL 
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APPENDIX  C 
PROGRAM  FOR  THE  3-D  GRAPH  USING  SURMAT  ABOUT  X-Y  PLANE 

MKICKKXXXKMMKKMKXKXIOOCXX   INFORMATION   xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

«  w 

*  THE  OBJECTIVE  OF  THIS  PROGRAM  IS  TO  CALCULATE  AND  DRAH  » 
«  THE  DIRECTIVITY  PATTERN  IN  THREE -DIMENSION.  « 

*  THE  DRAWING  METHOD  IS  THE  SUBROLTTINE  SURMAT  OF  DISSPLA.  * 
w  « 
xxKKXXXXXXXXXXXXxxxxxx  VARIABLE  DEFINITION  xxxxxxxxxxxxxxxxxxxxxxxxxxx 
w  * 
»  Z  :  DIRECTIVITY  IN  dB  » 
«  ZMM  :  DOUBLE  MATRIX  TO  GET  X-Y  PLANE  « 

*  2M  :  DUMMY  VARIABLE  OF  Z  » 
»  PI  :  PHASE  (DEGREE)  OF  SOURCE  A  » 
w  P2  :  PHASE  (DEGREE)  OF  SOURCE  B  » 

*  P3  :  PHASE  (DEGREE)  OF  SOURCE  C  » 
»  P4  :  VARIABLE  TO  CHANGE  DEGREE  TO  RADIAN  » 
«  P5  :  RADIAL  VALUE  OF  PI  » 
«  P6  :  RADIAL  VALUE  OF  P2  « 
«  P7  :  RADIAL  VALUE  OF  P3  » 
»  K  :  WAVE  NUMBER  » 

*  Dl  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  1  « 

*  D2  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  A  « 
«  D2  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  B     * 

*  D4  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  C  * 
»  PHI  :  PHI  IN  DEGREE  » 
»  I  :  THETA  IN  DEGREE  « 

*  PH  :  PHI  IN  RADIAN  « 
»  TH  :  THETA  IN  RADIAN  » 
»  PI  :  3.141S917  » 
»  » 
xxKKKXxxxxKXXXKKXxxxKX   VARIABLE  DECLARATION   mkxkxxxxxxkkxxxxxxkkxkxmkx 

DIMENSION  Z(  8100 ),ZMM( 90,90) 

REAL  P,T,R,X,Y,ZM,a,b,c 

REAL  P^,P5,P6,P7,W1,W2,W5,W^,W,G,H,K,D1,D2,D3,D^,PH,TH,PI 

INTEGER  PHI,P1,P2,P5,I 
c 
c 

PROGRAM  FOR  THE  UPPER  HEMISPHERE 

c 
C 

c — -~ — — — - — - — - — — 

C  A  PROGRAM  FOR  PLOTTING  BY  SHERPAOR  BY  TEK6I8 

C 

c 

CALL  COMPRS 
C     CALL  TEK618:  MOVE  C  TO  THE  UPPER  LINE  WHEN  DRAWING  WITH  TEK618 

CALL  RESET  ( 3HALL  ) 

CALL  PAGE(15. ,9.  ) 

CALL  PHYS0R(1.,3. ) 

CALL  SCMPLX 

CALL  AREA2D(1<*.,8.  ) 

CALL  X3NAME(1H  ,1) 

CALL  Y3NAME(1H  ,1) 

CALL  Z3NAME(1H  ,1) 

CALL  V0LM3D  ( 17. ,17. ,17. ) 

CALL  VUANGL(+20.,+20. ,500. ) 

CALL  INTAXS 

CALL  ZAXANG  (-90.  ) 

CALL  GRAF3D(0.,10.,80.,0.,10.,80.,0.,10.,40.  ) 

CALL  BLSUR 

CALL  BGNMAT( 90,90) 
C 
C 

C  INPUT  PARAMETERS 

C 

c 

PI =3.1415917 

K=PI 

Dl=l. 
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D2=l. 
D5=l. 
D<t=l. 
Pl=  180 


P2= 
P3= 
A  = 
B  = 
C  = 


0 

180 

1. 

0. 

0. 


c 

c 

C  MAIN  PROGRAM 

C 

c 

P^=PI/180 
P5=Pl*P<f 
P6=P2*P^ 
P7=P3*P^ 
DO  200  I=-90,0 
TH=I*P<t 

DO   100    PHI=1,360,4 
PH=PHI*P<t 

W1=K*D1*SIN( PH )»SIN( TH ) 
W2=K*D2*SIN( PH )*SIN(  TH )-P5 
H3=K*D3»SIN( TH )*COS(  PH )-P6 
W<+=K*D'+*SIN(  TH  )*COS(  PH  )+P7 
W=COS(  Wl  )+A»COS(  W2  )+B*COS(  W3  )+C*COS(  W<t ) 
G=SIN(  Nl  )-A*SIN(  W2  )-B*SIN(  W2  )+C*SIN(  W<t ) 
H=SQRT( W»W+G*G )/( 1+A+B+C ) 
21PHI)=20*ALOG10(H) 
IF(2(PHl).LE.-'+0.  )   THEN 

Z(PHI)=0. 
ELSE 

2(PHI)=2(PHI) 
END   IF 

R=<»0  . +Z(  PHI  ) 
P=  FLOAT(PHI) 
T=FLOAT(I) 

X=R*COS( TH  )»COS( PH )+40 . 
Y=R*COS(  TH  )»SIN(  PH  )+<^0. 
ZM=-R*SIN(TH) 
WRITE(8,1)  X,Y,ZM 
1  FORMAT(2X,3(E10.3,5X)) 

CALL  GETMAT(X,Y,2M,1,0) 
100      CONTINUE 
200    COfiTINUE 
C 
C 

C          A  PROGRAM  FOR  PLOTTING  BY  SHERPAOR  BY  TEK618 
C 

C 

CALL  ENDMAT(2MM,0) 

CALL  SURMAT  (ZMM, 1,90,1, 90, 0 ) 

CALL  DASH 

CALL  MARKER!  17) 

CALL  COMPLX 

CALL  HEIGHTt .2) 

CALL  SURVISl  'BOTTOM'  ) 

CALL  MESSAG( 'DIRECTIVITY  PATTERNS' ,19, 2.0, -.5) 

CALL  RESETt  'COMPLX' ) 

CALL  RESETC  'HEIGHT' ) 
c 
c 

PROGRAM  FOR  THE  LONER  HEMISPHERE 

c 
C 

CALL  BGNMATC 90,90) 
C 

C  INPUT  PARAMETERS 

C 

c 

PI=3.141S917 

K=PI 

01=1. 

D2=l. 


60 


D5=l. 
D^=l. 
Pl=  180 
P2=  0 
P3=  180 
A  =  1. 
B  =  0. 
C  =  0. 


C  MAIN  PROGRAM 

C 

c 

P^=PI/180 
PS=Pl«P<t 
P6=P2*P<t 
P7=P3*P<» 

DO  250  I=-180,-90 
TH=I*P^ 

DO  150  PHI=l,360,<<t 
PH=PHI*P^ 

W1=K»D1*SIN( PH )*COS( TH ) 
W2=K»02*SIN( PH  )*COS( TH  )*PS 
W3=K*D3*SIN( TH  )*COS( PH  )-P6 
H4=K*D4*SIN( TH )*COS( PH )+P7 
W=COS( Wl  )+A*COS( W2  )+B*COS( W3 )+C*COS( N4 ) 
G=SIN(  Wl  )-A*SIN(  W2  J+B*SIN(  W5  )-C*SIN(  W<t ) 
H=SQRT( W*H+G*G )/( 1+A+B+C ) 
Z(PHI)=20*AL0G10(HJ 
IF(Z(PHI).LE,-40.  )  THEN 

2(PHI)=0. 
ELSE 

Z(PHI)=Z(PHI) 
END  IF 

R=<*0.+Z(PHI) 
P=  FLOAT(PHI) 
T=FLOAT(I) 

X=R*COS( TH  )»COS(  PH )*40 . 
Y=R*COS( TH  )*SIN( PH )*40 . 
ZM=R*SIN(TH) 
WRITE(8,1)  X,Y,ZM 
1  FORMAT(2X,3(E10.3,SX)) 

CALL  GETMAT(X,Y,2M,1,0J 
ISO      CONTINUE 
250    CONTINUE 

CALL  ENDMAT(ZMM,0) 
CALL  SURMAT  (ZMM, 1,90,1, 90, 0 ) 
CALL  DASH 
CALL  MARKER(17) 
CALL  COMPLX 
CALL  HEIGHT( .2) 
CALL  SURVIS( 'BOTTOM' ) 
CALL  ENDPL(O) 
CALL  DONEPL 
CHECKSUM=012362 
STOP 
END 
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APPENDIX  D 
PROGRAM  FOR  THE  3-DIMENSIONAL  GRAPH  USING  VECTOR  DRAWING  MET 

MKXXKKMKKKXXXKKXXKKKKXM   INFORMATION   XXXXKXXXXXXXXKXXKXXXXXKMXXKKXXXXXK 

n  n 

»                    THE  OBJECTIVE  OF  THIS  PROGRAM  IS  TO  CALCULATE  AND  DRAW  « 

«        THE  DIRECTIVITY  PATTERN  IN  THREE-DIMENSION.  • 
»        THE  DRAWING  METHOD  IS  THE  VECTOR  DRAWING  METHOD  OF  DISSPLA.    * 

«  w 
KXKXKXXXXXXxxxxxxxxxxx  VARIABLE  DEFINITION   xxxxxxxxxxxxxxxxxxxxxxxxxxx 

«        Z  :  DIRECTIVITY  IN  dB  « 

«        XT  :  "TO"  POINT  OF  X  IN  DRAWING  VECTOR  » 

»        YT  :  "TO"  POINT  OF  Y  IN  DRAWING  VECTOR  » 

«        2T  :  "TO"  POINT  OF  Z  IN  DRAWING  VECTOR  « 

»        XF  :  "FROM"  POINT  OF  X  IN  DRAWING  VECTOR  « 

«        YF  :  "FROM"  POINT  OF  Y  IN  DRAWING  VECTOR  « 

«        ZF  :  "FROM"  POINT  OF  Z  IN  DRAWING  VECTOR  « 

«        PI  :  PHASE  (DEGREE)  OF  SOURCE  A  « 

«        P2  :  PHASE  (DEGREE)  OF  SOURCE  B  « 

«        P5  :  PHASE  (DEGREE)  OF  SOURCE  C  « 

»        P<*  :  VARIABLE  TO  CHANGE  DEGREE  TO  RADIAN  » 

«        P5  :  RADIAL  VALUE  OF  PI  « 

»        P6  :  RADIAL  VALUE  OF  P2  « 

»        P7  :  RADIAL  VALUE  OF  PZ  • 

«        K  :  WAVE  NUMBER  » 

«        01  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  1     » 

»        D2  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  A     » 

»        D2  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  B     » 

»        D<»  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  C      » 

»        PHI  :  PHI  IN  DEGREE  « 

»                    I  :  THETA  IN  DEGREE  « 

«        PH  :  PHI  IN  RADIAN  « 

»        TH  :  THETA  IN  RADIAN  » 

«        PI  :  3.141S917  » 

«  « 

KKMKKXKXKXKMXXWXKKKKXK   VARIABLE  DECLARATION   KKXXKKKKKMXKKKXXKXKKKXXKMX 

REAL  R,Z,XT,YT,XF,YF,ZF,ZT,a,b,c 

REAL  P<+,PS,P6,P7,Wl,W2,W3,W<t,W,G,H,K,Dl,D2,D3,D<»,PH,TH,PI 
INTEGER  PHI, PI, P2, PS, I 
C _ 

DIMENSION  Z(  0 :  2500  )  ,ZMM(  i*6  ,^6  )  ,V(  2000  )  ,U(  2000  ) 

REAL  P,T,R,X,Y,ZM,ZZM,A,B,C,V,U 

REAL  P'f,PS,P6,P7,Wl,W2,W5,W^,W,G,H,K,Dl,D2,D3,D4,PH,TH,PI 

INTEGER  PHI,P1,P2,P3,I,S 
C 
C 

C  INPUT  PARAMETERS 

C 


c 


PI=5. 1415917 

K=10.5S 

Dl=0.5 

D2=0.S 

D3=0.5 

04=0.5 

Pl=  0 

P2=180 

P3=180 

A  =  1. 

B  =  1. 

C  =  1. 


C                          MAIN  PROGRAM 
C 

c 

P4=PI/180 

P5=P1*P<+ 

P6=P2*P4 

P7=P3«P<f 

S=l 
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c 

CI 

c 

CI 
100 
200 
C 


DO  200  1=0,90 
TH=I*P^ 

DO  100  PHI =0,180 
PH=PHI»P4 
W1=K*D1*SIN( PH )»SIN( TH ) 

.   H2=K*D2*SIN(PH)*SIN(TH)-PS 

-   W5=K*D3*SIN(TH)*C0S(PH)-P6 
W4=K*D4*SIN( TH  )*COS( PH )+P7 
W=COS( Wl J+A»COS( H2  )+B*COS( W5 )+C»COS( W^ ) 
G=SIN(  HI  1-A*SIN(  W2  )-B*SIN(  W5  )+C*SIN(  W<* ) 
H=SQRT( W*W+G»G )/( 1+A+B+C  J 
2(PHI1=20*ALOG1Q(H) 
IF((2(PHI)  .GT.-IO.  ).AND.(2(PHI).LT.-9. ))  THEN 

Z(PHI)=  -20. 
ELSE 

2(PHI)=  2JPHI) 
END  IF 
2M=Z(PHI) 

WRITE(8,1)  I,PHI,ZM 
IF  (2M  .EQ.  -50.  )  THEN 
V(S)  =  FLOAT  (PHI) 
U(S)  =  FLOAT  (I) 
HRITE(8,1)  V(S),U(S),2M 
S=S+1 
END  IF 

WRITE! 8,1)  I, PHI ,2 
FORMAT(2X,2(I3,5X),F10.S) 

HRITE(8,1)  AXrL),BY(L),2M 
FORMATC  2X,3( F10.3,5X) ) 
CONTINUE 
CONTINUE 


A  PROGRAM  FOR  PLOTTING  BY  SHERPA  OR  BY  TEK618 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
STOP 
END 


,10) 
,10) 


MEDBUF 

COMPRS 

TEK618 

PHYSORd.  ,1.  ) 

PAGE( 15.5,12. ) 

AREA2D(10.5,9.1) 

FRAME 

COMPLX 

HEIGHTt .2) 

XNAMEt  'PHI(DEG)   $' 

YNAME( 'THETA(DEG)$' 

YTICKSC  2  ) 

XTICKS(  <¥ ) 

GRAF( 0,20,180,0,10,90) 

CURVE(V,U,S,-1) 

MESSAGJ  •THETA=2S  DEGREE ' ,15,11. ,8.5) 

MESSAG( 'K=S.5',5,11.,8.2) 

MESSAG( '01=1.  D2=l. ',11,11.,7.9) 

MESSAG( '03=1.  04=1. ' ,11, 11., 7. 6) 

MESSAG( 'P2=180  P3=0  P^=180 ' ,18,11. ,7.3 ) 

MESSAG(*A=1.  B=l.  C=l. ' ,14,11. ,7.0) 

RESET(  'HEIGHT'  ) 

HEIGHT!  .<¥) 

MESSAGt  ' TRACK ( -20. DB)' 

RESET( 'ALL'  ) 

NOCHEK 

END6R(0) 

PHYSORt 11.85,7.6) 

AREA20(3.5,2.5) 

FRAME 

ENDPL(O) 


,13,3.1,9.8) 
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APPENDIX  E 
PSD 


INFORMATION       XKMXXKXJOOOOOOtKKMXXKKKXKXKXKXKXXKK 


KXKXXXXKKXXKXXXXXXXKXXX 

» 

»  THE  OBJECTIVE  OF  THIS  PROGRAM  IS  TO  CALCULATE  AND  DRAW 

«  THE  DIRECTIVITY  PATTERN  IN  THREE -DIMENSION. 

»  THE  DRAWING  METHOD  IS  THE  SUBROLTTINE  SURMAT  OF  DISSPLA 

«  ABOUT  D-Q)  PLANE. 

M 

XXKXKKXKXXKXKKKXKKXXKX   VARIABLE  DEFINITION   XXXXXXKKXXXMKXKKKXK 

« 

DIRECTIVITY  IN  dB       _  _ 

:  DOUBLE  MATRIX  TO  GET  b-<P  plane 

DUMMY  VARIABLE  OF  THETA 

DUMMY  VARIABLE  OF  PHI 
;  DUMMY  VARIABLE  OF  DIRECTIVITY 
;  PHASE  (DEGREE)  OF  SOURCE  A 
:  PHASE  (DEGREE)  OF  SOURCE  B 
I  PHASE  (DEGREE)  OF  SOURCE  C 
:  VARIABLE  TO  CHANGE  DEGREE  TO  RADIAN 
;  RADIAL  VALUE  OF  PI 
;  RADIAL  VALUE  OF  P2 
;  RADIAL  VALUE  OF  P5 

WAVE  NUMBER 
:  DISTANCE  FROM  THE 
;  DISTANCE  FROM  THE 
;  DISTANCE  FROM  THE 
;  DISTANCE  FROM  THE 

:  PHI  IN  DEGREE 

THETA  IN  DEGREE 
;  PHI  IN  RADIAN 
1  THETA  IN  RADIAN 
:  2.1<tlS917 


« 
« 

KKKKKKXK 


« 

Z  : 

» 

2MM 

« 

X  : 

» 

Y  : 

« 

ZM 

» 

PI 

« 

P2 

« 

P5 

« 

P<¥ 

» 

P5 

« 

P6 

If 

P7 

« 

K  : 

» 

Dl 

« 

D2 

» 

D3 

» 

D** 

« 

PHI 

« 

I  : 

» 

PH 

« 

TH 

m 

PI  ; 

w 

GEOMETRICAL 
GEOMETRICAL 
GEOMETRICAL 
GEOMETRICAL 


CENTER  TO  THE  SOURCE 
CENTER  TO  THE  SOURCE 
CENTER  TO  THE  SOURCE 
CENTER  TO  THE  SOURCE 


« 

« 

« 
1  » 
A  » 
B  « 
C     » 

XXXXXXKX 


VARIABLE  DECLARATION   xxxxxxxxxxxxxxxxxx 
DIMENSION  Z(2S00),ZMM(<i6,<f6) 

REAL  X,Y,ZM,P<^>P5,P6,P7,Wl,W2,W3>W4,W,G,H,K,01,D2,D3,D<f,PH,TH,PI 
real  a>b>c 
INTEGER  PHI,PI,P2,P3,I 


A  PROGRAM  FOR  PLOTTING  BY  SHERPAOR  BY  TEK618 


CALL  COMPRS 

CALL  TEK618:  MOVE  C  TO  THE  UPPER  LINE  WHEN  DRAWING  WITH  TEK618 

CALL  RESET  (3HALL) 

CALL  PHYSOR(O.S,1. ) 

CALL  PAGE( 15.5,12.  ) 

CALL  SCMPLX 

CALL  AREA2D(10.5,9.15) 

CALL  X3NAME(  •THETA( DEGREE) ',13) 

CALL  Y3NAME( 'PHK DEGREE)' ,11) 

CALL  Z3NAME(  'DIRECTIVITY( DB  ) ' ,15  ) 

CALL  V0LM3D  ( 17. ,17. ,17.  ) 

CALL  VUANGL(+20.,+20.,S00.  ) 

CALL  INTAXS 

CALL  ZAXANG  ( -90. ) 

CALL  GRAF3D(0.,5.,90.,180.,6.,360.,-<*0.,2.,0.  ) 

CALL  BLSUR 

CALL  BGNMAT(<»6,46) 


INPUT  PARAMETERS 


PI  =3. 1<*  15917 

K=5.S 

Dl=l. 

D2=l. 

D3=l. 

D<»=1. 


64 


Pl=  180 
P2=  0 
P3=  180 
A  =  1. 
B  =  1. 
C  =  1. 


MAIN  PROGRAM 


100 
200 

C 


P^=PI/180 
PS=P1»P<» 
P6=P2*P4 
P7=P3*P4 
DO  200  1=0, 90 > 2 
TH=I«P<t 

DO  100  PHI=180,360,4 
PH=PHI*P<t 

W1=K*D1*SIN( PH )*SIN(  TH ) 
H2=K*D2*SIN( PH  )*SIN( TH )-PB 
H3=K*D2*SIN( TH  )*COS(  PH )-P6 
W<f=K*D4*SIN(  TH  )*COS(  PH  )+P7 
W=COS( Wl )+A*COS( W2  )+B*CQS( W3  )4C«C0S( W4 ) 
G=SIN( Wl  )-A«SIN( W2  )-B*SIN( W3 )+C»SIN( W4 ) 
H=SQRT( W*W+G*G  )/( 1+A+B+C ) 
Z(PHI)=20*ALaG10(H) 
IF(2(PHI).LE.-40. )  THEN 

Z(PHI)=-40. 
ELSE 

Z(PHIJ=  2(PHIJ 
END  IF 
X=I 
Y=PHI 
ZM=Z( PHI ) 
WRITE(8,1)  X,Y,2M 
F0RMAT(2X,3(E10.3,5X)) 
CALL  6ETMAT(X,Y,ZM,1,0) 
CONTINUE 
CONTINUE 


C 

C 

C 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 
CHECKSUM=012362 

STOP 

END 


A  PROGRAM  FOR  PLOTTING  BY  SHERPAOR  BY  TEK618 


ENDMAT(2MM,0) 

SURMAT  (ZMM, 1,^6,1, '♦6,0) 

DASH 

MARKER! 17) 

COMPLX 

2) 

BOTTOM ■ ) 

K=5.5' ,5,11. ,8.5) 

Dl=l.  D2=l. ",11,11. ,8.2) 

D3  =  l.  D<^=1.  ',11,11.  ,7.9) 

VIEW  ANGLE=  20,  20 ' ,18,11. ,7.6 ) 

P2  =  180  P3=0  P<f=180' ,18,11.,7.3) 

MESSAG( "Asl.  B=l.  C=l. ' ,14,11. ,7.0 ) 

RESETJ  'HEIGHT'  ) 

HEIGHT( .4) 

MESSAGl 'DIRECTIVITY  PATTERN' ,19,3.05,9.8) 

ENDGR(O) 

PHYSORI 11.35,7.6) 

AREA2D(3.5,2.5) 

FRAME 

DONEPL 


HEIGHT! 
SURVIS( 
MESSAG( 
MESSAGJ 
MESSAG( 
MESSAG( 
MESSAG( 
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APPENDIX  F 
TRACE 


MXKKMXXKKXKXKXKXXXXXXXK  INFORMATION  MXXXXXXXKXKKXXXXKKXXXXKKKXXXXXXXKK 
H  It 

»  THE  OBJECTIVE  OF  THIS  PROGRAM  IS  TO  CALCULATE  AND  DRAW        » 

«        THE  CONSTANT  dB  LEVEL  OF  THE  THREE-DIMENSIONAL  PLOT.  « 

»  THE  DRAWING  METHOD  IS  THE  SUBROUTINE  CURVE  OF  DISSPLA.  » 
»  » 

KKXXXXXXXXxxxxxxxxxxxx   VARIABLE  DEFINITION  xxxxxxxxxxxxxxxxxxxxxxxxxxx 

»  2  :  DIRECTIVITY  IN  dB  « 

«  V  :  ARRAY  TO  TRACK  JP  « 

»  u  :  ARRAY  TO  TRACK  V  K 

«  PI  :  PHASE  (DEGREE)  OF  SOURCE  A  » 

»  P2  :  PHASE  (DEGREE)  OF  SOURCE  B  « 

«  P3  :  PHASE  (DEGREE)  OF  SOURCE  C  * 

»  P<*  :  VARIABLE  TO  CHANGE  DEGREE  TO  RADIAN  » 

«  PS  :  RADIAL  VALUE  OF  PI  » 

*  P6  :  RADIAL  VALUE  OF  P2  » 
»  P7  :  RADIAL  VALUE  OF  P3  » 
»  K  :  WAVE  NUMBER  * 
»  Dl  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  1  » 
»  D2  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  A  * 

•  D3  :  DISTANCE"  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  B  « 
»  D4  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  C  » 
«  PHI  :  PHI  IN  DEGREE  » 
»  I  :  THETA  IN  DEGREE  * 
«  PH  :  PHI  IN  RADIAN  • 
»  TH  :  THETA  IN  RADIAN  « 
«  PI  :  3.1'H5917  « 

KKKKKKXKKKXXXXXXKXKKXK   VARIABLE  DECLARATION   MXXXMKKKMKMXKXWKMXKXKMKKKK 

DIMENSION  Z(  0 :  2500  )  >ZMM(  <*6 ,46  ) ,  V(  2000  )  ,U(  2000  ) 

REAL  P,T,R,X,Y,ZM,ZZM,A,B,C,V,U 

REAL  P'+,P5,P6,P7,W1,W2,W3,W4,W,G,H,K,D1,02,D5,D4,PH,TH,PI 

INTEGER  PHI,P1,P2,P3,I,S 
e 
C 

C  INPLTT  PARAMETERS 

C 

c 

PI=3,l<tlS917 

K=10.55 

D1=0.S 

D2=0,S 

03=0.5 

04=0.5 

Pl=    0 

P2=180 

P3=180 

A   =    1. 

B   =   1. 

C   =   1. 


C  MAIN  PROGRAM 

C 

c 

P4=PI/180 
P5=P1»P4 
P6=P2*P4 
P7=P3*P4 
S=l 

DO  200  1=0,90 
TH=I*P4 
DO  100  PHI=0,180 

PH=PHI»P4 

W1=K*D1*SIN( PH )*SIN( TH ) 

W2=K*D2*SIN( PH )«SIN( TH  )-P5 

W3=K*D3*SIN( TH  )*COS( PH )-P6 

W4=K*D4XSIN( TH )*COS( PH )+P7 

W=COS( Wl )+A*COS( W2 )+B*COS( W3 )+C»COS( W4 ) 
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c 

CI 

c 

CI 
100 
200 
C 


6=SIN( Wl )-A«SIN( H2 )-B«SIN( H3 )*C»SIN( H4 ) 

H=SQRT( W*W+G»G )/( 1+A+B+C ) 

Z(PHI)=20»AL0G10(H) 

1F((Z(PHI)  .GT.-IO, ).AND.(ZCPHI).LT.-9.))  THEN 

2(PHI)=  -50. 
ELSE 

Z(PHI)=  Z(PHI) 
END  IF 
ZM=Z(PHI) 

NRITEI8,1)  I,PHI,ZM 
IF  (2M  .EQ.  -30.  )  THEN 
V(S)  =  FLOAT  (PHI) 
U(S)  =  FLOAT  (I) 
WRITE(8,1)  V(S)>U(S)>ZM 
S=S+1 
END  IF 

HRITE(8,1)  I,PHI,Z 
FORMAT(2X,2(I2,SX),F10.5) 

WRITE(8,1)  AX(L),BY(L),ZM 
FORMAT(2X,5(F10.3,5X)) 
CONTINUE 
CONTINUE 


A  PROGRAM  FOR  PLOTTING  BY  SHERPA  OR  BY  TEK618 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
STOP 
END 


>10) 
,10) 


MEDBUF 

COMPRSjMOVE  C  TO  TEK619  TO  PRINT  USING  SHERPA 

TEK618 

PHYS0R(1.,1. ) 

PAGE( 15.5,12. ) 

AREA2Dtl0.5j9.1) 

FRAME 

COMPLX 

HEIGHT ( .2) 

XNAME(  'PHKDEG)   $' 

YNAME( 'THETAIDEG)*' 

YTICKSt  2 ) 

XTICKS( ^ ) 

GRAF(0,20,180,0,10,90) 

CURVE (V,U,S,-1) 

MESSAGt  'THETA=2S  DEGREE ' ,15,11. ,8.5) 

MESSAG( ■K=5.5' ,5,11. ,8.2) 

MESSAGI '01=1.  D2=l. ' ,11,11.,7.9) 

MESSAGt  'D3=l.  D'+sl.  '  ,11,11.  ,7.6  ) 

MESSAGI  'P2 =180  P3=0  P«^=180  '  ,18,11.  ,7.3  ) 

MESSAGt  'A=l.  B=l.  C=l. ' ,1^,11. ,7.0 ) 

RESET( 'HEIGHT' ) 

HEIGHT! .4) 

MESSAGt  ■TRACK(-20, 

RESETt  'ALL'  ) 

NOCHEK 

ENDGRtO) 

PHYSORt 11.85,7.6) 

AREA2Dt3.5,2.5) 

FRAME 

ENDPL(O) 


DB)', 13, 3. 1,9. 8) 
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APPENDIX  G 
LIST  OF  COMPUTER  OUTPUT 


•   SAMPLE  CASE   1 


ICXKKKKKXXXXXKKXKKKXKKKKMKXKKKXKKKKXMMKKKMKKKKMKKKKKKKKXXKKKKKKMXKKKXItXXKKMKKXXKKXICXXKKMKXKKlCXKK 

A  B  C  K  D,  D-,  D,  D,  (D-,      (D,         (p.         <P        6  H  F 

KXKKXKKXKKWXXKKXKKKKXKKKKKKXKXKKXKKXKXKKXKTtKXKMX>W(KKkKKKXK>fkKKMXKKXXXKXXXXXXKKKXKKKKXKKXXXKKXWM 


1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

31 

0.955567 

-0.394776 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

32 

0.970322 

-0.261679 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

33 

0.982069 

-0.157164 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

3't 

0.990826 

-0.080049 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

35 

0.996641 

-0.029226 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

36 

0.999579 

-0.003656 

1.0 

1.0 

1.0 

10.55 
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